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' Abstract The consistent quantum theory of self-phase modulation (SPM) and 

[ cross-phase modulation (XPM) for ultrashort light pulses (USPs) in 

■ medium with electronic Kerr-nonlinearity are developed. The approach 
{"S |' makes use of momentum operator of electrical field which takes account 

of the inertial behaviour of the nonlinearity. The spectrum of quantum 
C~| ■ fluctuations of squeezed-quadrature component as a function of response 

' time of nonlinearity and values of nonlinear phase shifts due to the SPM 

^ . and XPM effects, is investigated. 
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1. INTRODUCTION 

There are some methods to obtain USPs in a nonclassical state. One 
is the degenerate three-frequency parametric amplification which has 
high phase sensitivity. Other is SPM in a Kerr-nonlinear medium which 
in comparison with parametrical processes, does not require phase- 
matching. This property is a real advantage for the pulsed nonclas- 
sical state formation. The both participation of SPM and dispersion in 
the nonlinear medium leads to the optical soliton formation [1]. The 
quantum theory of SPM and XPM for USPs developed during last 15 
years have difficulties produced, in particular, by the multi-frequency 
structure of USPs and nonlinear kind of their interaction. 

The aim of our study is to develop the consistent quantum theory of 
nonlinear propagation of USPs. An account is taken of the role of a 
response time of the nonlinearity. The first attempt has been made in 
[2] to give quantum analysis of SPM of USPs for the case of relaxation 
nonlinearity of medium, and has been noted that the correct quantum 
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theory should take account of additional noise sources connected with 
nonlinear absorption. However, the theory developed in [2] does not take 
account of noise sources, so that the commutation relation for annihi- 
lation and creation photon operators is not fulfilled. This was carried 
out in the developed in [3] quantum theory of SPM, where noise sources 
are considered as a fluctuation addition to the relaxation nonlinearity 
of medium. Therefore, the results of [3] are connected with the range 
of carrier frequencies of the pulse, where the nonlinear absorption is im- 
portant. Actually the authors [3] have developed the quantum theory 
of USP propagation in a Raman active medium. In [3] the electronic 
Kerr nonlinearity with finite response time is modeled as the Raman 
active medium. Note, if we deal with the USPs propagation, for ex- 
ample, through fused-silica fibres only about 0.2 of the Kerr effect is 
attributable to the Raman oscillators and about 0.8 of the Kerr effect is 
due to electronic motion [4]. The case when the USP carrier frequency 
is far enough from some resonances (one- and two-photon and the Ra- 
man resonances), and therefore absorption is absent, was investigated by 
us in [5] where the finite response time is considered in the interaction 
Hamiltonian, and the commutation relation is exactly fulfilled. 

In the present work the results of the consistent quantum theory of 
SPM and XPM of USPs based on the use of the momentum operator 
(quantity of movement) are presented. The developed approach can be 
used when the response time is much shorter than the pulse duration and 
dispersion in nonlinear medium is neglected. There are no restrictions 
on the pulse intensity in our theory. 



2. QUANTUM THEORY OF SPM OF USPs 



The traditional way to describe the SPM of USP is based on the 
interaction Hamiltonian use, at that usually we solve the time-evolution 
equation. The transition to the spatial-evolution equation is realized 
substituting t with z/u, where z is the distance passed in medium and 
u is the group velocity. This approach seems to be good enough for 
single-mode radiation. If we deal with nonlinear propagation of USP 
then both t and z are present in analytical description. Thus, we use 
the momentum operator of pulse field related to the space-evolution [6] . 

We describe the SPM of USP using the momentum operator (cf. [5]) 



G sp m(z) = -hp / dt / H(t - h) N [n(t, z)h(h,z)] dt u (1.1) 
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where h(t,z) = A + (t, z)A(t, z), A + (t, z) (A(t,z)) is the Bose operator 
creating (annihilating) photons in a given cross-section z of the medium 
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at a given time t, Ti is the Planck's constant and N is the operator of 
normal ordering. The coefficient (3 is defined by the Kerr nonlinearity of 
a medium at stationary conditions. H{t) is the function of the nonlinear 
response of a medium; H{t) ^ at t > and H{t) = at t < 0. 
The expression under the first integral in (1.1) can be interpreted as 
generalized force acting in a defined cross-section z, which at the moment 
of time t depends only on the previous ones (t\ < t), i.e. it satisfies the 
causality principle. The nonlinear response function of a medium should 
be introduced as : 

H{t) = (l/rJe-V* (t > 0), (1.2) 

where r r is the response time of nonlinear medium, that we assume to 
be much shorter than the pulse duration r p . This nonlinearity takes 
place in absence of the one- and two- photons absorption and Raman 
resonances [1]. Besides, as it will be shown below, also in this limit case 
the account of finite nonlinear response time plays an important role. 
The space evolution equation for A(t,z) follows from (see [6]) 



..dA(t,z) 
-ih- 



dz 

and making use of (1.1), we finally get 



A(t,z),G spm {z) , (1.3) 



- i l -[3q [h(t, z)\ A(t, z) = 0, (1.4) 

where 

roo 

q[h{t,z)}= H(t 1 )[n(t-t 1 ,z)+n{t + t 1 ,z)]dt 1 . (1.5) 
Jo 

Here h(t, z) = A + (t, z)A(t, z) is the photon number "density". The Eq. 
(1.4) describes SPM in the moving coordinate frame: z = z , t = t —z/u, 
where t is the running time. According to (1.4) the operator n(t, z) does 
not depend on z, i.e. h(t,z) = h(t,z = 0) = ho(t) = Aq (t)Ao(t) (z = 
corresponds to the entrance into the medium). Taking account of this 
fact q[h(t,z)] = q[no(t)], and (1.5) can be rewritten as: 



q[n (t)]= fc(ti)fio(*-fi)ctti, (h(t) = H(\t\)). (1.6) 

Taking account of (1.6), we get the solution of (1.4) 

A(t, z) = exp {i iq [n (t)]}A (t). (1.7) 
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For hermitian conjugate operator of A(t, z) we have 

A+(t, z) = A+(t) exp {-i iq [n Q (t)}}. (1-8) 

In expressions (1.7), (1.8) 7 = (3z/2. For time-independent operator no, 
Eqs. (1.7), (1.8) give us the results for single mode radiation [7]. In the 
case of the instantaneous response of the nonlinearity H(t) = 5(t) we 
get [2] 

A(t,z) = e i2 ^ o{t) A (t), A + (t,z) = i f (t)e- i27 "' o(t) . (1.9) 

At the input to the medium Bose operators satisfy the commutation 
relation L4o(ii), Aq fe)] = 5(h —£2)- in a correct quantum methodology 
the analogical commutation relation must also be satisfied in nonlinear 
medium, i.e. for any z 

[A(t 1 ,z),A + (t 2 ,z)]=5(t 1 -t 2 ). (1.10) 

The calculation of average values of (1.9) and the normal ordering are 
followed by the non-integrable singularities appearance (for instance, 
5(t) appears in the exponent). These shortcomings are absent using 
(1.7) and (1.8). Besides, their use requires knowledge of an algebra of 
time-dependent Bose operators. The later has been developed in [2, 8]. 
The following permutation relations [8] 

i (*i)e d( * 2) = e^+^^-^Aoih), (1.11) 
e°^A+(h) = A+it^e 6 ^^^-^, (1.12) 

(where 0(t) = i^q[no(t)]) and the theorem of normal ordering [9] 



e o(*)=NexpW [e^W-1 



h (t-T r 6)de\, (1.13) 



are valid, where h(6) = e"^, 9 = t/r r . As [6(h), 6(t 2 )] = 0, then 
e 6(ti) e 6(t 2 ) = e 6(ti)+6(t 2 )_ Ma k in g use Q f re l a ti ns (1.11) and (1.12), 
one can verify that the commutation relation (1.10) is exactly fulfilled. 

As we analyse the SPM of an initial coherent USP, the equation on 
eigenvalues Ao(t)\a(t)) = a(t)\a(t)) is satisfied by Ao(t) (see [1]), where 
\ce(t)) is the initial coherent state, a(t) is the eigenvalue of Ao(t), and 
|a(t)| 2 = n-o(t) is the average photon density. Particular attention is 
given to the quantum fluctuations behaviour of X-quadrature; X(t, z) = 
[A(t, z) + A + (t, z)]/2. The behaviour of another quadrature component 
is shifted in phase with ir/2. For correlation function (see [5]) 

R(t, t) = (X(t, z)X(t + r, z)) - (X(t, z))(X(t + r, z)), (1.14) 
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where the brackets denote averaging over the initial coherent state of 
the pulse, making use of (1.11)-(1.13) we get 

R(t, r) = - [<5(t) - il>{t)h(r) sin 2$(t) + ^ 2 (t)g{r) sin 2 $(*)] . (1.15) 

In (1.15) the following notations have been introduced: $(t) = tp(t) + 
if(t), ip(t) = arg a(t) - the phase of USP or the heterodyne pulse's 
phase at balanced homodyne detection, ip(t) = 2~/no(t) = ij)Qp 2 (t) - 
the nonlinear phase shift resulted from SPM, ipo = tjj(0) = 27710 - the 
maximum nonlinear phase shift, p(t) - the envelope of USP (|ct(t)| = 
a p(t), p(0) = 1), and g(r) = J^h^hiO + r) old. The derivation 
of (1.15) took into consideration the fact that, in many experimental 
situations nonlinear phase shift per photon 7 < 1. The instantaneous 
spectral density of the quantum fluctuations of A-quadrature component 
Sx{w,t) = f^° 00 Rx(t,T)e lulT dT, according to (1.15), takes the form 

S x (n, t) = -[l- 2tp(t)L(n) sin 2$(t) + 4^ 2 (t)L 2 (Q) sin 2 *(t)] , (1.16) 

where L(O) = 1/[1 + O 2 ], O = ujt t . From (1.16) it follows that the 
quantum fluctuations level bellow the short noise one depends on the 
nonlinear phase shift ip{t) and &(t). At the initial phase of the pulse 
chosen optimal for a frequency = ujqt t - 
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/ \ 1 
ifo (t) = — arctan 



(i-i7) 



the spectral density (1.16) is 

s x (n , t) = 1{ [1 + ^ 2 (t)L 2 (n )] 1/2 - ii(t)L(n )} 2 . (lis) 

From (1.18) it follows that the spectral density adiabatically changes 
itself with the changing pulse's envelope and it is lower than 1/4 which 
corresponds to the coherent state of the pulse. From (1.18) we can also 
see that at the optimal phase ipo(t) with the increasing of the nonlinear 
phase shift ip(t), the spectral density Sx(&o,t) monotonously decreases. 
We point out that at r r = from (1.16) the results for single mode 
radiation can be obtained, at that the quantum fluctuation spectrum's 
level does not depend on the frequency. According to (1.18), increasing 
r r (increasing f^o) at fixed frequency Uo, the quadrature fluctuations 
level increases. 

3. QUANTUM THEORY OF XPM OF USPs 

We analyse now two-pulse propagation in an inertial nonlinear medium. 
We will consider that pulses have orthogonal polarizations or/and differ- 
ent frequencies. Then, besides SPM of each pulse, the XPM effect takes 
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place. Here we assume that the parametrical interaction of pulses can 
be neglected. In this case, the analysed process can be depicted making 
use of the following momentum operator [10]. 

G(z) = ]T GgL(z) + Gg$(z). (1.19) 

Here G&p m (z) is due to the SPM of j-pulse (j = 1,2) (see 1.1) and 
operator G X pm\Z) is connected with the XPM of pulses 

/oo rt 
dt / H(t-t 1 )[n 1 (t,z)n 2 (t 1 ,z)+n 2 (t,z)fi 1 (t 1 ,z)]dt 1 . 
-oo J —oo 

(1.20) 

In (1.20) hj(t,z) is the photon number operator of j-pulse and coef- 
ficients /3j, (3 are responsible for the SPM and XPM respectively. As 
stated earlier, the operator hj(t, z) commutes with G(z). In consequence 
hj(t, z) = hj(t, z = 0) = hoj(t), which means that the photons statistics 
of each pulse remains unchanged in the nonlinear medium. In accor- 
dance with (1.3), the evolution of A±(t, z) for the first pulse is given by 
(cf. (1.4)) 

dA f z ,Z) ~ i [\Piq [noM) + fa N, 2 (i)]} ii(t, z) = 0. (1.21) 

One can get for Ai(t, z) the similar equation changing the index 1 i — > 2 
in (1.21). The functions q[no,j{t)] is similar to (1.6). We remark that 
(1.4), (1.21), are written in the moving frame. The solution of (1.21) is 

ii(i,z) = e i Ti9["o,i(t)]+i7?[«o, 2 (t)]^ 0jl ^) ) (! 22) 

where 71 = [3\z/2, 7 = (3z. We define the correlation function of the in- 
vestigated light pulse as (1.14). For the spectra of quantum fluctuations 
of X\- quadrature we get 

S Xl (tl,t) = i{l-2^i(t)L(fi)sin2*i, 2 (t) 

+ 4[^(t)+Mt)Mt)]L 2 (n)sm 2 ^ 2 (t)}. (1.23) 

Here we denoted: $i, 2 (t) = + $i(*) = ipi(t) + (pi(t), rpi(t) = 

271^0,1 (*)) V'i (*) = 27noj(t) (j = 1, 2). As already mentioned, the phase 
$i(t) is connected with own parameters of the investigated pulse, and 
the phase ipi(t) with the XPM. From comparison of (1.23) and (1.16) 
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one can see that the XPM adds new terms in the multiplier and in the 
phase in (1.23). This circumstance allows us to control the fluctuations 
spectrum of the investigated pulse. At the initial phase 
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arctan 



Y>i(i) 



M 2 (t)L(n ) 



-Mt)-Mt), 



(1.24) 



chosen optimal for frequency Qq = ujQT r , the spectral density (1.23) is 

s Xl {n ,t) = i{i - 2L(n ) [^(t)+< 2 (t)L 2 (fio)] 1/2 -^ )2 (*)i(Oo)}, 

( L25 ) 

where ipl 2 (t) = ipi (t) + V'lWV^W- The influence of the second pulse to 
the spectrum of the first one is depicted in Figure 1. Increasing the pho- 




Figure 1 The spectrum of quantum fluctuations of the investigated light pulse at 
SPM and XPM as a function of reduced frequency Q. = uir r for different relations 
between average photon numbers of pulses no,2/fio,i: ( a )> 3 (b), 5 (c), 8 (d). The 
graphics are depicted for the case 71 = 72 = 2j, and nonlinear phase shift ipo,i(0) = 2 
at t = 0. At SPM only (curve (a)), the maximum squeezing is realized for SI = 0. 



ton number (intensity) of the second pulse, their quantum fluctuations 
can substantially increase the level of this ones for the investigated pulse. 
Changing the intensity of the second pulse, one can control the spectrum 
of quantum fluctuations of the the investigated pulse's quadrature. 

4. CONCLUSIONS 

The main result of our work is represented by the development of the 
consistent quantum theory of SPM and XPM of USPs in non-absorption 
nonlinear media, when the carrier frequencies of the pulse is far enough 
from any resonances. However, the consideration of the finite response 
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time of nonlinearity is necessary to get correct solutions ((1.7), (1.8), 
(1.22)). Our results are valid for response time of nonlinearity much 
shorter than duration of pulses but for unspecified intensities of USPs. 
In some sense, our results are complementary with the ones presented in 
[3] which treats the situation when the Raman resonance is important. 
The quantum theory of XPM of USPs in the developed approach here 
is presented for the first time. We have shown that the form and the 
level of the fluctuation spectrum can be controlled with the change of 
the pulse's phase and the another pulse's intensity in the presence of 
XPM. This fact can be used for quantum non-demolition measurements 
of parameters of the pulse [11]. 
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